home *** CD-ROM | disk | FTP | other *** search
/ CD ROM Paradise Collection 4 / CD ROM Paradise Collection 4 1995 Nov.iso / hobby / ast44src.zip / PLACALC2.C < prev    next >
C/C++ Source or Header  |  1995-02-11  |  33KB  |  859 lines
  1. /*
  2. ** Astrolog (Version 4.40) File: placalc2.c
  3. **
  4. ** IMPORTANT NOTICE: The graphics database and chart display routines
  5. ** used in this program are Copyright (C) 1991-1995 by Walter D. Pullen
  6. ** (astara@u.washington.edu). Permission is granted to freely use and
  7. ** distribute these routines provided one doesn't sell, restrict, or
  8. ** profit from them in any way. Modification is allowed provided these
  9. ** notices remain with any altered or edited versions of the program.
  10. **
  11. ** The main planetary calculation routines used in this program have
  12. ** been Copyrighted and the core of this program is basically a
  13. ** conversion to C of the routines created by James Neely as listed in
  14. ** Michael Erlewine's 'Manual of Computer Programming for Astrologers',
  15. ** available from Matrix Software. The copyright gives us permission to
  16. ** use the routines for personal use but not to sell them or profit from
  17. ** them in any way.
  18. **
  19. ** The PostScript code within the core graphics routines are programmed
  20. ** and Copyright (C) 1992-1993 by Brian D. Willoughby
  21. ** (brianw@sounds.wa.com). Conditions are identical to those above.
  22. **
  23. ** The extended accurate ephemeris databases and formulas are from the
  24. ** calculation routines in the program "Placalc" and are programmed and
  25. ** Copyright (C) 1989,1991,1993 by Astrodienst AG and Alois Treindl
  26. ** (alois@azur.ch). The use of that source code is subject to
  27. ** regulations made by Astrodienst Zurich, and the code is not in the
  28. ** public domain. This copyright notice must not be changed or removed
  29. ** by any user of this program.
  30. **
  31. ** Initial programming 8/28,30, 9/10,13,16,20,23, 10/3,6,7, 11/7,10,21/1991.
  32. ** X Window graphics initially programmed 10/23-29/1991.
  33. ** PostScript graphics initially programmed 11/29-30/1992.
  34. ** Last code change made 1/29/1995.
  35. */
  36.  
  37. #include "placalc.h"
  38.  
  39.  
  40. #ifdef PLACALC
  41. /*
  42. ** ---------------------------------------------------------------
  43. ** | Copyright Astrodienst AG and Alois Treindl, 1989,1991,1993  |
  44. ** | The use of this source code is subject to regulations made  |
  45. ** | by Astrodienst Zurich. The code is NOT in the public domain.|
  46. ** |                                                             |
  47. ** | This copyright notice must not be changed or removed        |
  48. ** | by any user of this program.                                |
  49. ** ---------------------------------------------------------------
  50. **
  51. ** Important changes:
  52. ** 11-jun-93 revision 1.12: fixed error which affected Mercury between -2100
  53. ** and -3100 (it jumped wildly).
  54. */
  55.  
  56. /***********************************************************
  57. ** $Header$
  58. **
  59. ** definition module for planetary elements
  60. ** and disturbation coefficients
  61. ** version HP-UX C  for new version with stored outer planets
  62. ** 31-jul-88
  63. ** by Alois Treindl 
  64. **
  65. ** ---------------------------------------------------------------
  66. ** | Copyright Astrodienst Zurich AG and Alois Treindl, 1989.    |
  67. ** | The use of this source code is subject to regulations made |
  68. ** | by Astrodienst Zurich. The code is NOT in the public domain.|
  69. ** |                                                             |
  70. ** | This copyright notice must not be changed or removed        |
  71. ** | by any user of this program.                                |
  72. ** ---------------------------------------------------------------
  73. **
  74. ***********************************************************/
  75.  
  76. /*
  77. ** In the elements degrees were kept as the units for the constants. This
  78. ** requires conversion to radians, when the actual calculations are performed.
  79. ** This approach is not the most efficient, but safer for development.
  80. ** Constant conversion could be done by writing all degree constants with
  81. ** value * DEGTORAD
  82. */
  83.  
  84. #define TIDAL_26  TRUE  /* decide wheter to use new or old lunar tidal
  85. term; a consistent system of delta t must be 
  86. used */
  87. #define MOON_TEST_CORR FALSE  /* to include more lunar terms in longitude */
  88.  
  89. REAL8 ekld[4] = {23.452294, -46.845, -.0059, 0.00181};
  90. /* ecliptic with epoch1900, Ekd(0..3) in basic */
  91.  
  92. struct eledata pd[MARS + 1] = {
  93. {/*earth*/  1.00000023, 365.25636042, EPOCH1900,
  94. 99.696678,  .9856473354,  1.089,    0,  
  95. 101.220833, 6189.03,  1.63,   0.012,
  96. 0.01675104, -0.00004180,  -0.000000126,
  97. 0,    0,    0,    0,
  98. 0,    0,    0},
  99. /*
  100. ** note 29 June 88 by Alois: G.M.Clemence, Astronomical Journal
  101. ** vol.53,p. 178 (1948) gives a correction to the perihel motion
  102. ** of -4.78" T, giving 6184.25 for the linear Term above. We have
  103. ** not yet applied this correction. It has been used in APAE 22,4
  104. ** on the motion of mars and does make an official impression.
  105. */
  106. {/*moon*/ 0.0025955307, 27.321661,  EPOCH1900,
  107. # if ! TIDAL_26 
  108. /*
  109. ** values from Improved Lunar Ephemeris, corresponding to tidal
  110. ** term -22.44"/cy and  consistent with delta t ~ 29.949 T*T
  111. */
  112. 270.4341638,  13.176396526808121, -4.08,  0.0068,
  113. # endif
  114. # if TIDAL_26 
  115. /*
  116. ** new values from Morrison 1979, with tidal term -26"/cy as
  117. ** stated in A.E. 1986 onwards, consistent with delta t ~ 44.3 T*T 
  118. ** correction: -1.54" + 2.33" T - 1.78" T*T
  119. */
  120. 270.4337361,  13.176396544528099, -5.86,  0.0068,
  121. # endif
  122. 334.329556, 14648522.52,  -37.17,   -0.045,
  123. 0.054900489,  0,    0,
  124. 259.183275, -6962911.23,  7.48 ,    0.008,
  125. 5.145388889,  0,    0},
  126. {/*mercury*/  .3870986, 87.969252,  EPOCH1900,
  127. 178.179078, 4.0923770233, 1.084,    0,
  128. 75.89969722,  5599.76,  1.061,    0,
  129. 0.20561421, .00002046,   -.000000030,
  130. 47.145944444, 4266.75,  .626,   0,
  131. 7.0028805555, 6.699,    -.066},
  132. {/*venus*/  .72333162,  224.700726, EPOCH1900,
  133. 342.767053, 1.6021687039, 1.1148,   0,
  134. 130.16383333, 5068.93,  -3.515,   0,
  135. 0.00682069, -.00004774, .000000091,
  136. 75.7796472223,3239.46,  1.476,    0,
  137. 3.3936305555, 3.621,    .0035},
  138. {/*mars*/ 1.5236914620, 686.9296097,  EPOCH1900,
  139. /* These are the corrected elements by Ross */
  140. 293.74762778, .524071163814,  1.1184,   0,
  141. 334.21820278, 6626.73,  .4675,    -0.0043,
  142. 0.09331290, .000092064, -.000000077,
  143. 48.786441667, 2775.57,  -.005,    -0.0192,
  144. 1.85033333, -2.430,   .0454}
  145. };
  146.  
  147. /*
  148. ** mimimum and maximum distances computed over 1000 years with plamimax,
  149. ** required for relative distances rgeo, where the distance is given
  150. ** as 100 when a planet is closest and as 0 when farthest from earth.
  151. */
  152. REAL8 rmima[CALC_N][2] = {  
  153.   { 0.98296342,  1.01704665},
  154.   { 0.00238267,  0.00271861},
  155.   { 0.54900496,  1.45169607},
  156.   { 0.26411287,  1.73597885},
  157.   { 0.37289847,  2.67626927},
  158.   { 3.94877993,  6.45627627},
  159.   { 7.99362824, 11.09276636},
  160.   {17.28622633, 21.10714104},
  161.   {28.81374786, 31.33507284},
  162.   {28.67716748, 50.29208774},
  163.   { 0.00,  0.00259553}, /* nodes don't get a real value */
  164.   { 0.00,  0.00259553},
  165.   { 7.36277475, 19.86585062}};
  166.  
  167.  
  168. struct sdat _sd [SDNUM] = {
  169.   114.50, 585.17493,
  170.   109.856,  191.39977,
  171.   148.031,  30.34583,
  172.   284.716,  12.21794,
  173.   114.508,  585.17656,
  174.   -0.56,    359.99213,
  175.   148.03, 30.34743,
  176.   284.72, 12.2196,
  177.   248.07, 1494.726615,
  178.   359.44, 359.993595,
  179.   109.86, 191.402867,
  180.   148.02, 30.348930,
  181.   114.503,  585.173715,
  182.   359.444,  359.989285,
  183.   148.021,  30.344620,
  184.   284.716,  12.21669,
  185.   148.0315, 30.34906264,
  186.   284.7158, 12.22117085,
  187.   220.1695,   4.284931111,
  188.   291.8024,   2.184704167
  189. };
  190.  
  191. REAL8 sa[SDNUM];
  192.  
  193. /*
  194. ** delta long = lampl * COS (lphase - arg) in seconds of arc
  195. ** delta rad  = rampl * COS (rphase - arg) in ninth place of log
  196. ** arg = j * sa (k) + i * ma (this planet)
  197. ** ma = mean anomaly
  198. ** sa = mean anomaly of disturbing planet, where this
  199. ** is taken from the aproximate value in sa[]
  200. ** For the COS (phase - arg) it is good enough to compute
  201. ** with 32 bit reals, because ampl and phase have only 
  202. ** four to five significant digits. 
  203. ** While saving constant space, it is costing execution time due
  204. ** to float/double conversions.
  205. **
  206. ** In basic, all correction terms for sun, mercury, venus and mars
  207. ** were contained in one array K(0..142,0..6); Nk(N,0) contained
  208. ** the index of the first term of planet N and Nk(N,1) the number
  209. ** of terms for this planet. Here, we use a  0 in the first column
  210. ** kor.j to indicate the end of the table for a planet.
  211. ** K(*) was a basic INTEGER array, therefore the amplitudes and phases
  212. ** had to be expressed as
  213. ** K(i,2) = ampl. of longitude in 0.001 seconds of arc
  214. ** K(i,3) = phase of longitude in 0.01 degrees
  215. ** K(i,4) = ampl. of radius in 9th place of log
  216. ** K(i,5) = phase of radius in 0.01 degrees.
  217. ** Here we have converted the amplitude of long. to seconds of arc
  218. ** and the phases to degrees.
  219. */
  220.  
  221. struct kor ARR earthkor[] = {  /* 11-jul-88 all terms to 0.020" longitude */
  222.   /*  j     i       lampl   lphase  rampl   rphase  k */
  223.   -1, 1,  0.013,  243,  28, 335,  8,  /* mercury */
  224.   -1, 3,  0.015,  357,  18, 267,  8,
  225.   -1, 4,  0.023,  326,  5,  239,  8,
  226.   -1, 0,  0.075,  296.6,  94, 205.0,  0,  /* Venus */
  227.   -1, 1,  4.838,  299.10, 2359, 209.08, 0,  
  228.   -1, 2,  0.074,  207.9,  69, 348.5,  0,
  229.   -1, 3,  0.009,  249,  16, 330,  0,
  230.   -2, 1,  .116, 148.90, 160,  58.40,  0,  
  231.   -2, 2,  5.526,  148.31, 6842, 58.32,  0,  
  232.   -2, 3,  2.497,  315.94, 869,  226.70, 0,  
  233.   -2, 4,  0.044,  311.4,  52, 38.8, 0,
  234.   -3, 2,  0.013,  176,  21, 90, 0,
  235.   -3, 3,  .666, 177.71, 1045, 87.57,  0,  
  236.   -3, 4,  1.559,  -14.75, 1497, 255.25, 0,  
  237.   -3, 5,  1.024,  318.15, 194,  49.50,  0,  
  238.   -3, 6,  0.017,  315,  19, 43, 0,
  239.   -4, 4,  .210, 206.20, 376,  116.28, 0,  
  240.   -4, 5,  .144, 195.40, 196,  105.20, 0,  
  241.   -4, 6,  .152, -16.20, 94, 254.80, 0,  
  242.   -5, 5,  0.084,  235.6,  163,  145.4,  0,
  243.   -5, 6,  0.037,  221.8,  59, 132.2,  0,
  244.   -5, 7,  .123, 195.30, 141,  105.40, 0,  
  245.   -5, 8,  .154, -.40, 26, 270.00, 0,  
  246.   -6, 6,  0.038,  264.1,  80, 174.3,  0,
  247.   -6, 7,  0.014,  253,  25, 164,  0,
  248.   -6, 8,  0.01, 230,  14, 135,  0,
  249.   -6, 9,  0.014,  12, 12, 284,  0,
  250.   -7, 7,  0.020,  294,  42, 203.5,  0,
  251.   -7, 8,  0.006,  279,  12, 194,  0,
  252.   -8, 8,  0.011,  322,  24, 234,  0,
  253.   -8, 12, 0.042,  259.2,  44, 169.7,  0,
  254.   -8, 14, 0.032,  48.8, 33, 138.7,  0,
  255.   -9, 9,  0.006,  351,  13, 261,  0,
  256.   1,  -1, .273, 217.70, 150,  127.70, 1,  /* mars */
  257.   1,  0,  0.048,  260.3,  28, 347,  1,
  258.   2,  -3, 0.041,  346,  52, 255.4,  1,
  259.   2,  -2, 2.043,  343.89, 2057, 253.83, 1,
  260.   2,  -1, 1.770,  200.40, 151,  295.00, 1,
  261.   2,  0,  0.028,  148,  31, 234.3,  1,  
  262.   3,  -3, .129, 294.20, 168,  203.50, 1,
  263.   3,  -2, .425, -21.12, 215,  249.00, 1,
  264.   4,  -4, 0.034,  71, 49, 339.7,  1,
  265.   4,  -3, .500, 105.18, 478,  15.17,  1,
  266.   4,  -2, .585, -25.94, 105,  65.90,  1,
  267.   5,  -4, 0.085,  54.6, 107,  324.6,  1,
  268.   5,  -3, .204, 100.80, 89, 11.00,  1,
  269.   6,  -5, 0.020,  186,  30, 95.7, 1,
  270.   6,  -4, .154, 227.40, 139,  137.30, 1,
  271.   6,  -3, .101, 96.30,  27, 188.00, 1,
  272.   7,  -5, 0.049,  176.5,  60, 86.2, 1,
  273.   7,  -4, .106, 222.70, 38, 132.90, 1,
  274.   8,  -5, 0.052,  348.9,  45, 259.7,  1,
  275.   8,  -4, 0.021,  215.2,  8,  310,  1,
  276.   8,  -6, 0.010,  307,  15, 217,  1,
  277.   9,  -6, 0.028,  298,  34, 208.1,  1,
  278.   9,  -5, 0.062,  346,  17, 257,  1,
  279.   10, -6, 0.019,  111,  15, 23, 1,
  280.   11, -7, 0.017,  59, 20, 330,  1,
  281.   11, -6, 0.044,  105.9,  9,  21, 1,
  282.   13, -8, 0.013,  184,  15, 94, 1,
  283.   13, -7, 0.045,  227.8,  5,  143,  1,
  284.   15, -9, 0.021,  309,  22, 220,  1,
  285.   17, -9, 0.026,  113,  0,  0,  1,
  286.   1,  -2, .163, 198.60, 208,  112.00, 2,  /* jupiter */
  287.   1,  -1, 7.208,  179.53, 7067, 89.55,  2,
  288.   1,  0,  2.600,  263.22, 244,  -21.40, 2,
  289.   1,  1,  0.073,  276.3,  80, 6.5,  2,
  290.   2,  -3, 0.069,  80.8, 103,  350.5,  2,
  291.   2,  -2, 2.731,  87.15,  4026, -2.89,  2,
  292.   2,  -1, 1.610,  109.49, 1459, 19.47,  2,
  293.   2,  0,  0.073,  252.6,  8,  263,  2,
  294.   3,  -3, .164, 170.50, 281,  81.20,  2,
  295.   3,  -2, .556, 82.65,  803,  -7.44,  2,
  296.   3,  -1, .210, 98.50,  174,  8.60, 2,
  297.   4,  -4, 0.016,  259,  29, 170,  2,
  298.   4,  -3, 0.044,  168.2,  74, 79.9, 2,
  299.   4,  -2, 0.080,  77.7, 113,  347.7,  2,
  300.   4,  -1, 0.023,  93, 17, 3,  2,
  301.   5,  -2, 0.009,  71, 14, 343,  2,
  302.   1,  -2, 0.011,  105,  15, 11, 3,  /* saturn */
  303.   1,  -1, .419, 100.58, 429,  10.60,  3,
  304.   1,  0,  .320, 269.46, 8,  -7.00,  3,
  305.   2,  -2, .108, 290.60, 162,  200.60, 3,
  306.   2,  -1, .112, 293.60, 112,  203.10, 3,
  307.   3,  -2, 0.021,  289,  32, 200.1,  3,
  308.   3,  -1, 0.017,  291,  17, 201,  3,
  309.   ENDMARK
  310. };
  311.  
  312. struct kor ARR mercurykor[] = {
  313.   1,  -1, .711, 35.47,  491,  305.28, 4,
  314.   2,  -3, .552, 161.15, 712,  71.12,  4,
  315.   2,  -2, 2.100,  161.15, 2370, 71.19,  4,
  316.   2,  -1, 3.724,  160.69, 899,  70.49,  4,
  317.   2,  0,  .729, 159.76, 763,  250.00, 4,
  318.   3,  -3, .431, 105.37, 541,  15.53,  4,
  319.   3,  -2, 1.329,  104.78, 1157, 14.84,  4,
  320.   3,  -1, .539, 278.95, 14, 282.00, 4,
  321.   4,  -2, .484, 226.40, 234,  136.02, 4,
  322.   5,  -4, .685, -10.43, 849,  259.51, 4,
  323.   5,  -3, 2.810,  -10.14, 2954, 259.92, 4,
  324.   5,  -2, 7.356,  -12.22, 282,  255.43, 4,
  325.   5,  -1, 1.471,  -12.30, 1550, 77.75,  4,
  326.   5,  0,  .375, -12.29, 472,  77.70,  4,
  327.   2,  -1, .443, 218.48, 256,  128.36, 5,
  328.   4,  -2, .374, 151.81, 397,  61.63,  5,
  329.   4,  -1, .808, 145.93, 13, 35.00,  5,
  330.   1,  -1, .697, 181.07, 708,  91.38,  6,
  331.   1,  0,  .574, 236.72, 75, 265.40, 6,
  332.   2,  -2, .938, 36.98,  1185, 306.97, 6,
  333.   2,  -1, 3.275,  37.00,  3268, 306.99, 6,
  334.   2,  0,  .499, 31.91,  371,  126.90, 6,
  335.   3,  -1, .353, 25.84,  347,  295.76, 6,
  336.   2,  -1, .380, 239.87, 0,  0,  7,
  337.   ENDMARK
  338. };
  339.  
  340. struct kor ARR venuskor[] = {
  341.   -1, 2,  .264, -19.20, 175,  251.10, 8,
  342.   -2, 5,  .361, 167.68, 55, 77.20,  8,
  343.   1,  -1, 4.889,  119.11, 2246, 29.11,  9,
  344.   2,  -2, 11.261, 148.23, 9772, 58.21,  9,
  345.   3,  -3, 7.128,  -2.57,  8271, 267.42, 9,
  346.   3,  -2, 3.446,  135.91, 737,  47.37,  9,
  347.   4,  -4, 1.034,  26.54,  1426, 296.49, 9,
  348.   4,  -3, .677, 165.32, 445,  75.70,  9,
  349.   5,  -5, .330, 56.88,  510,  -33.36, 9,
  350.   5,  -4, 1.575,  193.93, 1572, 104.21, 9,
  351.   5,  -3, 1.439,  138.08, 162,  229.90, 9,
  352.   6,  -6, .143, 84.40,  236,  -5.80,  9,
  353.   6,  -5, .205, 44.20,  256,  314.20, 9,
  354.   6,  -4, .176, 164.30, 70, 75.70,  9,
  355.   8,  -5, .231, 180.00, 25, 75.00,  9,
  356.   3,  -2, .673, 221.62, 717,  131.60, 10,
  357.   3,  -1, 1.208,  237.57, 29, 149.00, 10,
  358.   1,  -1, 2.966,  208.09, 2991, 118.09, 11,
  359.   1,  0,  1.563,  268.31, 91, -7.60,  11,
  360.   2,  -2, .889, 145.16, 1335, 55.17,  11,
  361.   2,  -1, .480, 171.01, 464,  80.95,  11,
  362.   3,  -2, .169, 144.20, 250,  54.00,  11,
  363.   ENDMARK
  364. };
  365.  
  366. struct kor ARR marskor[] = {
  367.   -1, 1,  .115, 65.84,  684,  156.14, 12,
  368.   -1, 2,  .623, 246.03, 812,  155.77, 12,
  369.   -1, 3,  6.368,  57.60,  556,  -32.06, 12,
  370.   -1, 4,  .588, 57.24,  616,  147.28, 12,
  371.   -2, 5,  .138, 39.18,  157,  309.39, 12,
  372.   -2, 6,  .459, 217.58, 82, 128.10, 12,
  373.   -1, -1, .106, 33.60,  141,  303.45, 13,
  374.   -1, 0,  .873, 34.34,  1112, 304.05, 13,
  375.   -1, 1,  8.559,  35.10,  6947, 304.45, 13,
  376.   -1, 2,  13.966, 20.50,  2875, 113.20, 13,
  377.   -1, 3,  1.487,  22.18,  1619, 112.38, 13,
  378.   -1, 4,  .175, 22.46,  225,  112.15, 13,
  379.   -2, 2,  .150, 18.96,  484,  266.42, 13,
  380.   -2, 3,  7.355,  158.64, 6412, 68.62,  13,
  381.   -2, 4,  4.905,  154.09, 1985, 244.70, 13,
  382.   -2, 5,  .489, 154.39, 543,  244.50, 13,
  383.   -3, 3,  .216, 111.06, 389,  21.10,  13,
  384.   -3, 4,  .355, 110.64, 587,  19.17,  13,
  385.   -3, 5,  2.641,  280.58, 2038, 190.60, 13,
  386.   -3, 6,  .970, 276.06, 587,  6.75, 13,
  387.   -3, 7,  .100, 276.20, 116,  6.40, 13,
  388.   -4, 5,  .152, 232.48, 259,  142.60, 13,
  389.   -4, 6,  .264, 230.47, 387,  139.75, 13,
  390.   -4, 7,  1.156,  41.64,  749,  312.67, 13,
  391.   -4, 8,  .259, 37.92,  205,  128.80, 13,
  392.   -5, 8,  .172, -8.99,  234,  260.70, 13,
  393.   -5, 9,  .575, 164.48, 308,  74.60,  13,
  394.   -6, 10, .115, 113.70, 145,  23.53,  13,
  395.   -6, 11, .363, 285.69, 144,  196.00, 13,
  396.   -7, 13, .353, 48.83,  85, 319.10, 13,
  397.   -8, 15, 1.553,  170.14, 110,  81.00,  13,
  398.   -8, 16, .148, 170.74, 154,  259.94, 13,
  399.   -9, 17, .193, 293.70, 23, 22.80,  13,
  400.   1,  -3, .382, 46.48,  521,  316.25, 14,
  401.   1,  -2, 3.144,  46.78,  3894, 316.39, 14,
  402.   1,  -1, 25.384, 48.96,  23116,  318.87, 14,
  403.   1,  0,  3.732,  -17.62, 1525, 117.81, 14,
  404.   1,  1,  .474, -34.60, 531,  59.67,  14,
  405.   2,  -4, .265, 192.88, 396,  103.12, 14,
  406.   2,  -3, 2.108,  192.72, 3042, 102.89, 14,
  407.   2,  -2, 16.035, 191.90, 22144,  101.99, 14,
  408.   2,  -1, 21.869, 188.35, 16624,  98.33,  14,
  409.   2,  0,  1.461,  189.66, 1478, 279.04, 14,
  410.   2,  1,  .167, 191.04, 224,  280.81, 14,
  411.   3,  -4, .206, 167.11, 338,  76.13,  14,
  412.   3,  -3, 1.309,  168.27, 2141, 76.24,  14,
  413.   3,  -2, 2.607,  228.41, 3437, 139.74, 14,
  414.   3,  -1, 3.174,  207.20, 1915, 115.83, 14,
  415.   3,  0,  .232, 207.78, 240,  298.06, 14,
  416.   4,  -4, .178, 127.25, 322,  36.16,  14,
  417.   4,  -3, .241, 200.69, 389,  110.02, 14,
  418.   4,  -2, .330, 267.57, 413,  179.86, 14,
  419.   4,  -1, .416, 221.88, 184,  128.17, 14,
  420.   1,  -2, .155, -38.20, 191,  231.58, 15,
  421.   1,  -1, 1.351,  -34.10, 1345, 235.85, 15,
  422.   1,  0,  .884, 288.05, 111,  39.90,  15,
  423.   1,  1,  .132, 284.88, 144,  15.67,  15,
  424.   2,  -2, .620, 35.15,  869,  305.30, 15,
  425.   2,  -1, 1.768,  32.50,  1661, 302.51, 15,
  426.   2,  0,  .125, 18.73,  103,  119.90, 15,
  427.   3,  -2, .141, 47.59,  199,  318.06, 15,
  428.   3,  -1, .281, 40.95,  248,  310.75, 15,
  429.   ENDMARK
  430. };
  431.  
  432. struct m45dat m45[NUM_MOON_CORR] = {
  433.   /*    l,     l',  F,  D,      Long,    Lat,         Par),*/
  434.   { 0,   0, 0,   4,    13.902,     14.06,  0.2607},
  435.   { 0,   0, 0,   2,    2369.912,   2373.36,     28.2333},
  436.   { 1,   0, 0,   4,     1.979,      6.98,  0.0433},
  437.   { 1,   0, 0,   2,     191.953,    192.72,  3.0861},
  438.   { 1,   0, 0,   0,   22639.500,  22609.1,     186.5398},
  439.   { 1,   0, 0,  -2,   -4586.465,  -4578.13,     34.3117},
  440.   { 1,   0, 0,  -4, -38.428,    -38.64,  0.6008},
  441.   { 1,   0, 0,  -6,  -0.393,     -1.43,  0.0086},
  442.   { 0,   1, 0,   4,  -0.289,     -1.59, -0.0053},
  443.   { 0,   1, 0,   2, -24.420,    -25.10, -0.3000},
  444.   { 0,   1, 0,   0,    -668.146,   -126.98,     -0.3997},
  445.   { 0,   1, 0,  -2,    -165.145,   -165.06,      1.9178},
  446.   { 0,   1, 0,  -4,  -1.877,     -6.46,  0.0339},
  447.   { 0,   0, 0,   3,   0.403,     -4.01,  0.0023},
  448.   { 0,   0, 0,   1,    -125.154,   -112.79,     -0.9781},
  449.   { 2,   0, 0,   4,   0.213,      1.02,  0.0054},
  450.   { 2,   0, 0,   2,  14.387,     14.78,  0.2833},
  451.   { 2,   0, 0,   0, 769.016,    767.96, 10.1657},
  452.   { 2,   0, 0,  -2,    -211.656,   -152.53,     -0.3039},
  453.   { 2,   0, 0,  -4, -30.773,    -34.07,  0.3722},
  454.   { 2,   0, 0,  -6,  -0.570,     -1.40,  0.0109},
  455.   { 1,   1, 0,   2,  -2.921,    -11.75, -0.0484},
  456.   { 1,   1, 0,   0,    -109.673,   -115.18,     -0.9490},
  457.   { 1,   1, 0,  -2,    -205.962,   -182.36,      1.4437},
  458.   { 1,  1,  0,  -4,  -4.391,     -9.66,  0.0673},
  459.   { 1,  -1, 0,  4,    0.283,      1.53,  0.0060},
  460.   { 1,  -1, 0,  2,   14.577,     31.70,  0.2302},
  461.   { 1,  -1, 0,  0,  147.687,    138.76,  1.1528},
  462.   { 1,  -1, 0,  -2,  28.475,     23.59, -0.2257},
  463.   { 1,  -1, 0,  -4,   0.636,      2.27, -0.0102},
  464.   { 0,  2,  0,  2,   -0.189,     -1.68, -0.0028},
  465.   { 0,  2,  0,  0,   -7.486,     -0.66, -0.0086},
  466.   { 0,  2,  0,  -2,  -8.096,    -16.35,  0.0918},
  467.   { 0,  0,  2,  2,   -5.741,     -0.04, -0.0009},
  468.   { 0,  0,  2,  0,  -411.608,    -0.2,      -0.0124},
  469.   { 0,  0,  2,  -2, -55.173,    -52.14, -0.1052},
  470.   { 0,  0,  2,  -4,   0.025,     -1.67,  0.0031},
  471.   { 1,  0,  0,  1,   -8.466,    -13.51, -0.1093},
  472.   { 1,  0,  0,  -1,  18.609,      3.59,  0.0118},
  473.   { 1,  0,  0,  -3,   3.215,      5.44, -0.0386},
  474.   { 0,  1,  0,  1,   18.023,     17.93,  0.1494},
  475.   { 0,  1,  0,  -1,   0.560,      0.32, -0.0037},
  476.   { 3,  0,  0,  2,    1.060,      2.96,  0.0243},
  477.   { 3,  0,  0,  0,   36.124,     50.64,  0.6215},
  478.   { 3,  0,  0,  -2, -13.193,    -16.40, -0.1187},
  479.   { 3,  0,  0,  -4,  -1.187,     -0.74,  0.0074},
  480.   { 3,  0,  0,  -6,  -0.293,     -0.31,  0.0046},
  481.   { 2,  1,  0,  2,   -0.290,     -1.45, -0.0051},
  482.   { 2,  1,  0,  0,   -7.649,    -10.56, -0.1038},
  483.   { 2,  1,  0,  -2,  -8.627,     -7.59, -0.0192},
  484.   { 2,  1,  0,  -4,  -2.740,     -2.54,  0.0324},
  485.   { 2,  -1, 0,  2,    1.181,      3.32,  0.0213},
  486.   { 2,  -1, 0,  0,    9.703,     11.67,  0.1268},
  487.   { 2,  -1, 0,  -2,  -2.494,     -1.17, -0.0017},
  488.   { 2,  -1, 0,  -4,   0.360,      0.20, -0.0043},
  489.   { 1,  2,  0,  0,   -1.167,     -1.25, -0.0106},
  490.   { 1,  2,  0,  -2,  -7.412,     -6.12,  0.0484},
  491.   { 1,  2,  0,  -4,  -0.311,     -0.65,  0.0044},
  492.   { 1,  -2, 0,  2,    0.757,      1.82,  0.0112},
  493.   { 1,  -2, 0,  0,    2.580,      2.32,  0.0196},
  494.   { 1,  -2, 0,  -2,   2.533,      2.40, -0.0212},
  495.   { 0,  3,  0,  -2,  -0.344,     -0.57,  0.0036},
  496.   { 1,  0,  2,  2,   -0.992,     -0.02,  0},
  497.   { 1,  0,  2,  0,  -45.099,     -0.02, -0.0010},
  498.   { 1,  0,  2,  -2,  -0.179,     -9.52, -0.0833},
  499.   { 1,  0,  -2, 2,   -6.382,     -3.37, -0.0481},
  500.   { 1,  0,  -2, 0,   39.528,     85.13, -0.7136},
  501.   { 1,  0,  -2, -2,   9.366,      0.71,     -0.0112},
  502.   { 0,  1,  2,  0,    0.415,      0.10,    0.0013},
  503.   { 0,  1,  2,  -2,  -2.152,     -2.26, -0.0066},
  504.   { 0,  1,  -2, 2,   -1.440,     -1.30,  0.0014},
  505.   { 0,  1,  -2, -2,   0.384,      0.0,   0.0},
  506.   { 2,  0,  0,  1,   -0.586,     -1.20, -0.0100},
  507.   { 2,  0,  0,  -1,   1.750,      2.01,  0.0155},
  508.   { 2,  0,  0,  -3,   1.225,      0.91,     -0.0088},
  509.   { 1,  1,  0,  1,    1.267,      1.52,  0.0164},
  510.   { 1,  -1, 0,  -1,  -1.089,      0.55,      0},
  511.   { 0,  0,  2,  -1,   0.584,      8.84,  0.0071},
  512.   { 4,  0,  0,  0,    1.938,      3.60,  0.0401},
  513.   { 4,  0,  0,  -2,  -0.952,     -1.58, -0.0130},
  514.   { 3,  1,  0,  0,   -0.551,      0.94, -0.0097},
  515.   { 3,  1,  0,  -2,  -0.482,     -0.57, -0.0045},
  516.   { 3,  -1, 0,  0,    0.681,      0.96,      0.0115},
  517.   { 2,  0,  2,  0,   -3.996,      0,   0.0004},
  518.   { 2,  0,  2,  -2,   0.557,     -0.75,     -0.0090},
  519.   { 2,  0,  -2, 2,   -0.459,     -0.38, -0.0053},
  520.   { 2,  0,  -2, 0,   -1.298,      0.74,      0.0004},
  521.   { 2,  0,  -2, -2,   0.538,      1.14, -0.0141},
  522.   { 1,  1,  -2, -2,   0.426,      0.07,     -0.0006},
  523.   { 1,  -1, 2,  0,   -0.304,      0.03,  0.0003},
  524.   { 1,  -1, -2, 2,   -0.372,     -0.19, -0.0027},
  525.   { 0,  0,  4,  0,    0.418,      0,   0},
  526.   { 2,  -1, 0,  -1,  -0.352,     -0.37,     -0.0028}
  527. };
  528.  
  529. # if MOON_TEST_CORR
  530. /* moon additional correction terms */
  531. struct m5dat {    
  532.   REAL8 lng;
  533.   int  i0,i1,i2,i3;
  534. } m5[] = {
  535.   /*    lng,  l,     l',  F,  D,      */
  536.   0.127,  0,  0,  0,  6,
  537.   -0.151, 0,  2,  0,  -4,
  538.   -0.085, 0,  0,  2,  4,
  539.   0.150,  0,  1,  0,  3,
  540.   -0.091, 2,  1,  0,  -6,
  541.   -0.103, 0,  3,  0,  0,
  542.   -0.301, 1,  0,  2,  -4,
  543.   0.202,  1,  0,  -2, -4,
  544.   0.137,  1,  1,  0,  -1,
  545.   0.233,  1,  1,  0,  -3,
  546.   -0.122, 1,  -1, 0,  1,
  547.   -0.276, 1,  -1, 0,  -3,
  548.   0.255,  0,  0,  2,  1,
  549.   0.254,  0,  0,  2,  -3,
  550.   -0.100, 3,  1,  0,  -4,
  551.   -0.183, 3,  -1, 0,  -2,
  552.   -0.297, 2,  2,  0,  -2,
  553.   -0.161, 2,  2,  0,  -4,
  554.   0.197,  2,  -2, 0,  0,
  555.   0.254,  2,  -2, 0,  -2,
  556.   -0.250, 1,  3,  0,  -2,
  557.   -0.123, 2,  0,  2,  2,
  558.   0.173,  2,  0,  -2, -4,
  559.   0.263,  1,  1,  2,  0,
  560.   0.130,  3,  0,  0,  -1,
  561.   0.113,  5,  0,  0,  0,
  562.   0.092,  3,  0,  2,  -2,
  563.   0,  99, 0,  0,  0 /* end mark */
  564. };
  565. # endif /* MOON_TEST_CORR */
  566.  
  567.  
  568. /*****************************************************
  569. $Header: deltat.c,v 1.10 93/01/27 14:37:06 alois Exp $
  570. deltat.c
  571. deltat(t): returns delta t (in julian days) from universal time t
  572. is included by users
  573. ET = UT +  deltat
  574.  
  575. ---------------------------------------------------------------
  576. | Copyright Astrodienst Zurich AG and Alois Treindl, 1989.    |
  577. | The use of this source code is subject to regulations made  |
  578. | by Astrodienst Zurich. The code is NOT in the public domain.|
  579. |                                                             |
  580. | This copyright notice must not be changed or removed        |
  581. | by any user of this program.                                |
  582. ---------------------------------------------------------------
  583.  
  584. ******************************************************/
  585.  
  586. double deltat(jd_ad)
  587. double jd_ad;
  588. { 
  589.   static short int ARR dt[] = { /* in centiseconds */
  590.   /*
  591.   ** dt from 1637 to 2000, as tabulated in A.E. 
  592.   ** the values 1620 - 1636 are not taken, as they fit
  593.   ** badly the parabola 25.5 t*t for the next range. The
  594.   ** best crossing point to switch to the parabola is
  595.   ** 1637, where we have fitted the value for continuity
  596.   */
  597.   6780, 6500, 6300, 
  598.   6200, 6000, 5800, 5700, 5500,
  599.   5400, 5300, 5100, 5000, 4900,
  600.   4800, 4700, 4600, 4500, 4400,
  601.   4300, 4200, 4100, 4000, 3800, /* 1655 - 59 */
  602.   3700, 3600, 3500, 3400, 3300,
  603.   3200, 3100, 3000, 2800, 2700,
  604.   2600, 2500, 2400, 2300, 2200,
  605.   2100, 2000, 1900, 1800, 1700,
  606.   1600, 1500, 1400, 1400, 1300,
  607.   1200, 1200, 1100, 1100, 1000,
  608.   1000, 1000, 900,  900,  900,
  609.   900,  900,  900,  900,  900,
  610.   900,  900,  900,  900,  900,  /* 1700 - 1704 */
  611.   900,  900,  900,  1000, 1000,
  612.   1000, 1000, 1000, 1000, 1000,
  613.   1000, 1000, 1100, 1100, 1100,
  614.   1100, 1100, 1100, 1100, 1100,
  615.   1100, 1100, 1100, 1100, 1100,
  616.   1100, 1100, 1100, 1100, 1200, /* 1730 - 1734 */
  617.   1200, 1200, 1200, 1200, 1200,
  618.   1200, 1200, 1200, 1200, 1300,
  619.   1300, 1300, 1300, 1300, 1300,
  620.   1300, 1400, 1400, 1400, 1400,
  621.   1400, 1400, 1400, 1500, 1500,
  622.   1500, 1500, 1500, 1500, 1500, /* 1760 - 1764 */
  623.   1600, 1600, 1600, 1600, 1600,
  624.   1600, 1600, 1600, 1600, 1600,
  625.   1700, 1700, 1700, 1700, 1700,
  626.   1700, 1700, 1700, 1700, 1700,
  627.   1700, 1700, 1700, 1700, 1700,
  628.   1700, 1700, 1600, 1600, 1600, /* 1790 - 1794 */
  629.   1600, 1500, 1500, 1400, 1400,
  630.   1370, 1340, 1310, 1290, 1270, /* 1800 - 1804 */
  631.   1260, 1250, 1250, 1250, 1250,
  632.   1250, 1250, 1250, 1250, 1250,
  633.   1250, 1250, 1240, 1230, 1220,
  634.   1200, 1170, 1140, 1110, 1060,
  635.   1020, 960,  910,  860,  800,
  636.   750,  700,  660,  630,  600,  /* 1830 - 1834 */
  637.   580,  570,  560,  560,  560,
  638.   570,  580,  590,  610,  620,
  639.   630,  650,  660,  680,  690,
  640.   710,  720,  730,  740,  750,
  641.   760,  770,  770,  780,  780,
  642.   788,  782,  754,  697,  640,  /* 1860 - 1864 */
  643.   602,  541,  410,  292,  182,
  644.   161,  10, -102, -128, -269,
  645.   -324, -364, -454, -471, -511,
  646.   -540, -542, -520, -546, -546,
  647.   -579, -563, -564, -580, -566,
  648.   -587, -601, -619, -664, -644, /* 1890 - 1894 */
  649.   -647, -609, -576, -466, -374,
  650.   -272, -154, -2, 124,  264,
  651.   386,  537,  614,  775,  913,
  652.   1046, 1153, 1336, 1465, 1601,
  653.   1720, 1824, 1906, 2025, 2095,
  654.   2116, 2225, 2241, 2303, 2349, /* 1920 - 1924 */
  655.   2362, 2386, 2449, 2434, 2408,
  656.   2402, 2400, 2387, 2395, 2386,
  657.   2393, 2373, 2392, 2396, 2402,
  658.   2433, 2483, 2530, 2570, 2624,
  659.   2677, 2728, 2778, 2825, 2871,
  660.   2915, 2957, 2997, 3036, 3072, /* 1950 - 1954 */
  661.   3107, 3135, 3168, 3218, 3268,
  662.   3315, 3359, 3400, 3447, 3503,
  663.   3573, 3654, 3743, 3829, 3920,
  664.   4018, 4117, 4223, 4337, 4449,
  665.   4548, 4646, 4752, 4853, 4959,
  666.   5054, 5138, 5217, 5296, 5379, /* 1980 - 1984 */
  667.   5434, 5487, 5532, 5582, 5630, /* 1985 - 89 from AE 1991 */
  668.   5686, 5757, 5900, 5900, 6000, /* AE 1993 and extrapol */
  669.   6050, 6100, 6150, 6200, 6250, /* 1995 - 1999 */
  670.   6300};          /* 2000 */
  671.   double yr, cy, delta;
  672.   long iyr, i;
  673.   yr = (jd_ad + 18262) / 365.25 + 100.0;    /*  year  relative 1800 */
  674.   cy = yr / 100;
  675.   iyr =  (long) (RFloor(yr) + 1800);   /* truncated to integer, rel 0 */
  676. #if TIDAL_26    /* Stephenson formula only when 26" tidal
  677.   term in lunar motion */
  678.   if (iyr >= 1637  && iyr < 2000) {
  679.   i = iyr - 1637;
  680.   delta = dt[i] * 0.01 + (dt[i+1] - dt[i]) * (yr - RFloor(yr)) * 0.01;
  681.   } else if (iyr >= 2000) { /* parabola, fitted at value[2000] */
  682.   delta = 25.5 * cy * cy  - 25.5 * 4 + 63.00; 
  683.   } else if (iyr >= 948) {  /* from 948 - 1637 use parabola */
  684.   delta = 25.5 * cy * cy; 
  685.   } else {  /* before 984 use other parabola */ 
  686.   delta = 1361.7  + 320 * cy + 44.3 * cy * cy;  /* fits at 948 */
  687.   }
  688. #else    /* use Clemence value + 5 sec before 1690, new dt afterwards */
  689.   cy -= 1;  /* epoch 1900 */
  690.   if (iyr >= 1690  && iyr < 2000) {
  691.   i = iyr - 1637;
  692.   delta = dt[i] * 0.01 + (dt[i+1] - dt[i]) * (yr - RFloor(yr)) * 0.01;
  693.   } else if (iyr >= 2000) { /* parabola, fitted at value[2000] */
  694.   delta = 29.949 * cy * cy  - 29.949 * 4 + 63.0;  
  695.   } else {  
  696.   delta = 5 + 24.349 + 72.3165 * cy + 29.949 * cy * cy; /* fits at 1690 */
  697.   }
  698. #endif
  699.   return delta / 86400.0;
  700. }
  701.  
  702.  
  703. /*******************************************
  704. $Header: d2l.c,v 1.9 91/11/16 16:24:20 alois Exp $
  705.  
  706. double to long with rounding, no overflow check
  707. *************************************/ 
  708. long d2l(x)
  709. double x;
  710. {
  711.   if (x >=0)
  712.     return ((long) (x + 0.5));
  713.   else
  714.     return (-(long) (0.5 - x));
  715. }
  716.  
  717.  
  718. /*
  719. * $Header$
  720. *
  721. * A collection of useful functions for centisec calculations.
  722.  
  723. ---------------------------------------------------------------
  724. | Copyright Astrodienst Zurich AG and Alois Treindl, 1991.    |
  725. | The use of this source code is subject to regulations made  |
  726. | by Astrodienst Zurich. The code is NOT in the public domain.|
  727. |                                                             |
  728. | This copyright notice must not be changed or removed        |
  729. | by any user of this program.                                |
  730. ---------------------------------------------------------------
  731. *******************************************************/
  732.  
  733. double degnorm(p)
  734. double p;
  735. {
  736.   if (p < 0) 
  737.     do {
  738.       p += 360.0;
  739.     } while (p < 0);
  740.   else if (p >= 360.0)
  741.     do {
  742.       p -= 360.0;
  743.     } while (p >= 360.0);
  744.   return (p);
  745. }
  746.  
  747.  
  748. /*********************************************************
  749. $Header: julday.c,v 1.9 91/11/16 16:25:06 alois Exp $
  750. *********************************************************/
  751.  
  752. /*
  753. ** This function returns the absolute Julian day number (JD)
  754. ** for a given calendar date.
  755. ** The aruguments are a calendar date: day, month, year as integers,
  756. ** hour as double with decimal fraction.
  757. ** If gregflag = 1, Gregorian calendar is assumed, gregflag = 0
  758. ** Julian calendar is assumed.
  759. **
  760. ** The Julian day number is system of numbering all days continously
  761. ** within the time range of known human history. It should be familiar
  762. ** for every astrological or astronomical programmer. The time variable
  763. ** in astronomical theories is usually expressed in Julian days or
  764. ** Julian centuries (36525 days per century) relative to some start day;
  765. ** the start day is called 'the epoch'.
  766. ** The Julian day number is a double representing the number of
  767. ** days since JD = 0.0 on 1 Jan -4712, 12:00 noon.
  768. ** Midnight has always a JD with fraction .5, because traditionally
  769. ** the astronomical day started at noon.
  770. **
  771. ** NOTE: The Julian day number is named after the monk Julianus. It must
  772. ** not be confused with the Julian calendar system, which is named after
  773. ** Julius Cesar, the Roman politician who introduced this calendar.
  774. **
  775. ** Original author: Marc Pottenger, Los Angeles.
  776. ** with bug fix for year < -4711   15-aug-88 by Alois Treindl
  777. **
  778. ** References: Oliver Montenbruck, Grundlagen der Ephemeridenrechnung,
  779. **             Verlag Sterne und Weltraum (1987), p.49 ff
  780. **
  781. ** related functions: revjul() reverse Julian day number: compute the
  782. **             calendar date from a given JD
  783. */
  784.  
  785. double julday(month, day, year, hour, gregflag)
  786. int month;
  787. int day;
  788. int year;
  789. double hour;
  790. int gregflag;
  791. {
  792.   double jd, u, u0, u1, u2;
  793.  
  794.   u = year;
  795.   if (month < 3)
  796.     u -=1;
  797.   u0 = u + 4712.0;
  798.   u1 = month + 1.0;
  799.   if (u1 < 4)
  800.     u1 += 12.0;
  801.   jd = RFloor(u0*365.25)
  802.     + RFloor(30.6*u1+0.000001)
  803.     + day + hour/24.0 - 63.5;
  804.   if (gregflag) {
  805.     u2 = RFloor(ABS8(u) / 100) - RFloor(ABS8(u) / 400);
  806.     if (u < 0.0)
  807.       u2 = -u2;
  808.     jd = jd - u2 + 2;            
  809.     if ((u < 0.0) && (u/100 == RFloor(u/100)) && (u/400 != RFloor(u/400)))
  810.       jd -= 1;
  811.   }
  812.   return jd;
  813. }
  814.  
  815.  
  816. /*********************************************************
  817. $Header: revjul.c,v 1.9 91/11/16 16:25:37 alois Exp $
  818. *********************************************************/
  819.  
  820. /*
  821. ** revjul() is the inverse function to julday(), see the description there.
  822. ** Arguments are julian day number, calendar flag (0=julian, 1=gregorian)
  823. ** return values are the calendar day, month, year and the hour of
  824. ** the day with decimal fraction (0 .. 23.999999).
  825. **
  826. ** Original author Mark Pottenger, Los Angeles.
  827. ** with bug fix for year < -4711 16-aug-88 Alois Treindl
  828. */
  829.  
  830. void revjul(jd, gregflag, jmon, jday, jyear, jut)
  831. double jd;
  832. int gregflag;
  833. int *jmon;
  834. int *jday;
  835. int *jyear;
  836. double *jut;
  837. {
  838.   double u0, u1, u2, u3, u4;
  839.  
  840.   u0 = jd + 32082.5;
  841.   if (gregflag) {
  842.     u1 = u0 + RFloor(u0/36525.0) - RFloor(u0/146100.0) - 38.0;
  843.     if (jd >= 1830691.5) u1 +=1;
  844.       u0 = u0 + RFloor(u1/36525.0) - RFloor(u1/146100.0) - 38.0;
  845.   }
  846.   u2 = RFloor(u0 + 123.0);
  847.   u3 = RFloor((u2 - 122.2) / 365.25);
  848.   u4 = RFloor((u2 - RFloor(365.25 * u3)) / 30.6001);
  849.   *jmon = (int)(u4-1.0);
  850.   if (*jmon > 12)
  851.     *jmon -= 12;
  852.   *jday = (int)(u2 - RFloor(365.25 * u3) - RFloor(30.6001 * u4));
  853.   *jyear = (int)(u3 + RFloor((u4 - 2.0) / 12.0) - 4800);
  854.   *jut = (jd - RFloor(jd + 0.5) + 0.5) * 24.0;
  855. }
  856. #endif /* PLACALC */
  857.  
  858. /* placalc2.c */
  859.